Mid-infrared absorption in YBa 2 Cu 3 06: 
Failure of spin-wave theory in undoped cuprates? 
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The optical conductivity a{ui) of undoped YBa2Cus06 is studied in detail in the mid-infrared 
range. Substitutions on all but the Ba site are used to identify the prominent absorption processes 
at 2800 and 3800 cm" 1 . Experimental evidence for bimagnon-plus-phonon absorption is collected. 
A more critical analysis of the lineshape and the spectral weight reveals the limits of this approach. 
Although phonon-2-magnon multiple scattering seems to reproduce the lineshape, the necessary 
coupling is unrealistically large. The strong increase of high frequency spectral weight with increasing 
temperature makes the failure of spin- wave theory even more evident. 

PACS numbers: 74.72.-h 74.25.Gz, 78.30.Hv, 75.40. Gb, 75.50.Ee 



The undoped parent compounds of the high T c 
cuprates are regarded as an almost ideal realization of 
a two-dimensional (2D) spin 1/2 Heisenberg antiferro- 
magnet. Despite the low dimensionality and the low spin 
the excitations are thought to be spin-waves with a well- 
defined dispersion fjl], as opposed to e.g. ID systems, 
where a spinon continuum is observed in neutron scatter- 
ing Q). In the cuprates, a spin- wave dispersion has been 
extracted throughout the whole Brillouin zone from the 
maxima in neutron scattering intensities, but energies are 
rather high, large backgrounds are observed and the mag- 
nitude of quantum corrections is unclear Moreover, 
the assumption that magnons are not well-defined parti- 
cles at the Brillouin zone boundary was a keypoint in the 
successful description of the photoemission data of insu- 
lating Sr 2 Cu0 2 Cl 2 in Ref. {§. Two-magnon (2M) Ra- 
man scattering shows several anomalies in the cuprates, 
in particular a very broad lineshape, spectral weight at 
high energies and a finite signal in Ai g geometry ||. A 
large body of theoretical work has been dedicated to this 
problem, and the importance of resonance phenomena jf| 
was emphasized. Other treatments include the interac- 
tion with phonons |t]] and extensions of the Heisenberg 
model H . Certainly, the strong influence of the charge 
transfer (CT) resonance on the Raman spectra compli- 
cates the problem significantly. 

Optical spectroscopy probes the magnetic excitations 
more directly. The main peak in mid-infrared ab- 
sorption (MIR) spectra of La 2 Cu04 and other single 
layer cuprates has been interpreted by Lorenzana 
and Sawatzky [ [10| in terms of bimagnon-plus-phonon 
(BIMP) absorption. A similar feature was reported in 
the bilayer system YBa 2 Cu 3 6 (YBC0 6 ) Q, in 2D 
S=l La 2 Ni0 4 @@ and in ID S=l/2 Sr 2 Cu0 3 @@. 
Good agreement is achieved in ID |C3 because quan- 
tum fluctuations are included ab initio, and for the 2D 
S=l nickelates |l^] because fluctuations beyond spin- 
wave theory are small. However, in the cuprates the 
estimate for the spectral weight of BIMPs is one order 



of magnitude too small. Moreover, a large amount of 
spectral weight is observed above the BIMP absorption. 
Interpretations in terms of multi-magnon-plus-phonon 
absorption |M , d-d transitions Jl2| and CT excitons 
have been proposed. We challenge these approaches and 
suggest that a full account of our MIR data in the un- 
doped cuprates has to include quantum fluctuations be- 
yond spin-wave theory. This might provide an important 
feedback to the Raman experiment and the basic picture 
of the undoped cuprates. 

Specific ionic substitutions help to identify the promi- 
nent resonances in the MIR spectrum of YBC06- In par- 
ticular we use oxygen isotope substitution to distinguish 
vibrational from electronic degrees of freedom and to give 
direct experimental evidence for the existence of BIMP 
absorption. Substituting rare earth (RE) elements for Y 
increases the lattice parameter a and thereby changes the 
exchange constant J. The dependence of the BIMP fre- 
quency on J (a) is found to be similar to the one reported 
for 2M Raman scattering ]il||i~7| ]. 

Single crystals of YBa 2 Cii307 were grown in BaZrOa 
(BZO) crucibles. Crystals grown using this technique 
exhibit a superior purity (> 99.995 at. %) @|9). The 
full exchange of the oxygen isotope was described else- 
where [ po[ . For comparison, single crystals grown in 
Y 2 3 stabilized Zr0 2 (YSZ) crucibles f|| were stud- 
ied as well. A finite amount of Y in the single crys- 
tals of RE . 8 Yo. 2 Ba 2 Cu 3 6 (RE=Pr, Gd) studied here 
is due to the Y 2 3 in the Zr0 2 crucibles. In the case 
of YBa 2 Cu 3 _ y Zn !y 6 |§, we estimate 0.05 < y < 0.07 
pof . The oxygen content of all crystals was fixed to a 
value 06 by annealing them for about a week in a flow of 
high purity Argon (99.998 %) at 750°C or in ultra high 
vacuum at 700°C. We calculated a(oj) by inverting the 
Fresnel equations for the experimentally measured trans- 
mission and reflection data. Measurements were carried 
out with the electric field vector polarized parallel and 
perpendicular to the afo-plane. The small remnants of 
interference fringes in some of the calculated spectra of 
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FIG. 1. a(u) of YBC0 6 at 4 and 300 K. Dashed and dot- 
ted lines: In- and inter-plane contributions to BIMP absorp- 
tion for J = 780 cm -1 , J12/J = 0.1 and hujph = 530 cm -1 . 
Dash-dotted line: sum of the two. Left inset: 2M Raman and 
IR BIM frequencies as a function of Jw/J- Right inset: J 
and hujph as a function of J12/J. 



<t(uj) are artifacts caused by deviations of the measured 
data from the assumed ideal case of absolutely flat and 
plane parallel surfaces. 

We first collect evidence for the applicability of a 
magnon-phonon interpretation of the main resonance. In 
Fig. 1 we display cr(ui) of YBCOe up to the onset of 
CT absorption for T=4 and 300 K. Note the very low 
values of cr(u>), which are two (four) orders of magni- 
tude lower than for YBCOg.i (YBCO7) in this frequency 
range. Following the interpretation of the single layer 
compounds |l(J we ascribe the main peak at 2800 cm -1 
to BIMP absorption. A magnetic origin is also favored 
by recent measurements of the pressure dependence of 
the MIR and Raman spectrum of Sr2Cu02Cl2 [p3j . In 
the cuprates bimagnon absorption is forbidden due to 
inversion symmetry. It only becomes weakly allowed by 
symmetry breaking effects like impurities or the combi- 
nation with a phonon. Let us compare the experimental 
data with predictions of BIMP theory on a qualitative 
basis (Fig. 2). The two upper plots show spectra of two 
samples of YBCOg grown in different crucibles (top panel 
YSZ, below BZO). The strong temperature dependence 
of the BIMP peak frequency is most likely due to the re- 
duction of spin stiffness with increasing temperature (left 
inset, open/full symbols YSZ/BZO). The width 7 of the 
BIMP resonance is smaller in the cleaner sample grown 



in a BZO crucible (right inset) |24|. We therefore expect 
that impurity scattering is not negligible in the deter- 
mination of the correct lineshapes. Substituting Zn on 
Cu sites indeed broadens the 2800 cm -1 peak drastically 
(panel 3 from top in Fig. 2). A finite phonon contri- 
bution to the BIMP peak is evident from the frequency 
shift induced by oxygen isotope substitution in YBCOg 
(see inset of same panel). The measured isotope shift of 
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FIG. 2. a(u)) for different substitutions and polarizations. 
Diamonds and broken lines in the top panel are theoretical 
results. See text for details. Insets top panel: T dependence 
of peak frequency and width of the 2800 cm" 1 peak in YBCOe 
samples grown in YSZ (open circles, top panel) and BZO (full 
circles, second panel). Inset panel 2: Spectral weight of the 
2800 and 3800 cm" 1 peaks. 



28 ± 8 cm" 1 is consistent with the BIMP interpretation, 
assuming that the longitudinal stretching phonon of ap- 
proximately 550 - 600 cm" 1 is excited. Substitution of Y 
with Pr or Gd leads to a significant frequency shift of the 
main peak. Similar shifts were observed in 2M Raman 
scattering jl6],[l7) and were explained by the dependence 
of J on the lattice parameter a. Similar to the case of 
Zn, the disorder on the Y site after substitution of 80% 
of Pr or Gd enhances 7. Finally, the 2800 cm" 1 peak is 
not observed if the electric field is applied parallel to the 
c axis, in agreement with the single layer data ||. The 
spikes on top of the main resonance and the other sharp 
features have not yet been identified. A relation to BIMP 
absorption is unlikely, as their position is identical in the 
RE samples. 

For a more critical analysis of peak frequency, line- 
shape and spectral weight we extend the BIMP theory to 
the bilayer case. In order to obtain the coupling to light 
we start from a Heisenberg Hamiltonian which takes into 
account a dependence of the in- and inter-plane exchange 



2 



constants J and J12 on the external electric field E and 
the phonon coordinates fin]: 



H 



E 

L=l,2 <i,j> 



J2 J(E, ujSi.iSij + ^ J 12 (E, u)S M S 



'2,/ 



where z and j label nearest neighbor Cu sites in a 2D 
square lattice, L labels the two planes in a single bilayer, 
and u denotes the displacements of O ions. Only Einstein 
phonons are considered. The different phonons modulate 
the intersite hopping and the on-site energies on both 
Cu and O sites. Modulations are taken into account to 
second order since we also included phonon-2M multi- 
ple scattering processes in a refined approach (see be- 
low). In Ref. fnj a high energy (HE) approximation was 
used which is inappropriate for zone center excitations. 
Our RPA results for two interacting magnons show that 
the HE approximation reproduces the line shape rather 
well. Both RPA and HE produce 2M bound states with 
a strong dispersion in momentum space. However, the 
sharp resonances at the low energy side of the BIMP re- 
ported in Ref. Jh} are removed in RPA. 

Both the Raman and infrared 2M peak frequencies 
were calculated as a function of j = J12 / J (left inset of 
Fig. 1). In the infrared case, the phonon frequency hu} p h 
still has to be added. At T=4K the experimental BIMP 
and 2M Raman spectra peak at 2795 and 2720 ±10 cm" 1 
H , respectively. From these we can deduce the values of 
J and hujph for a given ratio J12/ J (right inset of Fig. 1). 
We obtain J = 790±10 cm" 1 and j = 0.08±0.04. Neutron 
scattering suggests j=0.1 - 0.15 |E5 26| and fiw p j, w 550 
- 600 cm" 1 [pTj for the relevant longitudinal stretching 
phonon mode. However, a finite next-nearest neighbor 
coupling J' will shift the values of J and J12 consider- 
ably @. 

Let us compare the calculated and measured lineshapes 
and oscillator strengths of the BIMP peak. The calcu- 
lated BIMP absorption for J=780 cm" 1 , Ji 2 =0.1J and 
huj p h = 530 cm" 1 is plotted together with the experimen- 
tal curve in Fig. 1 (dash-dotted line). An offset of 0.1 
51" 1 cm" 1 has been used. In a bilayer we have to distin- 
guish two contributions: a photon can flip two spins (i) in 
the same layer (dashed line, in-plane) or (ii) in adjacent 
layers (dotted line, inter-plane). For (i) a rough estimate 
of the 2M binding energy in the Ising limit is J. This im- 
plies that the 2M energy at the zone boundary is about 
J below the upper cut-off of the 2M-spectrum, which ap- 
proximately corresponds to the BIMP position in Fig. 1. 
Similarly, the binding energy is J12 for case (ii) which 
explains the maximum of the calculated inter-plane con- 
ductivity just below the 2M cut-off. The estimated rel- 
ative spectral weight of inter- and in-plane contributions 
is 0.06 - 0.3 for j—0.1. Due to this small value the spec- 
tral weight is similar in YBCOg and in La2Cu04 both 
experimentally and theoretically. The perturbatively es- 
timated spectral weight is a factor of 8 - 15 too small 
compared to experiment p9j. 





peak A 


peak B 


A/B 


2M cut-off 


YBa 2 Cu 3 06 


2800 


3800 


0.74 


3700 = 4.72J 


La2Cu04 


3300* 


4500* 


0.73 


4500 = 4.63J 


Sr 2 Cu0 2 Cl 2 


2900* 


4000* 


0.73 


4000 = 4.63J 



TABLE I. Measured frequencies of the two main MIR ab- 
sorption peaks, their ratio, and the 2M cut-off (all frequencies 
in cm -1 ). * taken from Ref. [9]. 



Regarding the lineshape the calculated curve is sharper 
than the experimental one, but the width of the BIMP 
peak is sample dependent, as stated above. As in the sin- 
gle layer cuprates, the real problem is obviously at higher 
frequencies: the strong peak at 3800 cm -1 remains un- 
explained. It is likely that the high energy anomaly has 
the same origin in MIR and Raman spectra. This is sub- 
stantiated by the absence of the anomaly in both spectro- 
scopies in S=l La2Ni04 p2] , pp| . In several cuprates, the 
frequency ratio of the two dominant MIR peaks is about 
0.73 (Table 1), strongly suggesting a common magnetic 
origin. In other terms, the second peak in both single and 
bilayers is close to the 2M cut-off. One way of shifting 
spectral weight to the 2M cut-off is to consider a finite in- 
teraction between phonons and bimagnons. A dimension- 
less coupling constant for phonon-2M multiple scattering 
is defined as X P 2M = jj (^r ) (u 2 ), from which we estimate 
X V 2M ~ —0.02... + 0.01 for the stretching phonon mode 
IPIIJ32I , A negative value of X P 2M translates into a re- 
pulsive phonon-2M interaction and shifts spectral weight 
to higher frequencies. A way to test the reliability of our 
estimate of X P 2M is to compare the linear coupling dJ/du 
with the experimental pressure dependence of J [p3| . 
There, our estimate is 1 - 2 times smaller. However, an 
excellent fit to the data (diamonds in top panel of Fig. 2) 
is obtained only if we assume X P 2M = —0.2 (dashed line) 
and add the BIMP contribution of the apical stretching 
phonon. The apical contribution is expected to have a 
5 times smaller weight and a negligible phonon- 2M cou- 
pling (dotted line). This large value of X P 2M makes such 
a scenario very unlikely. Phonon-magnon scattering pro- 
cesses have also been considered for the explanation of 
the width of the 2M Raman resonance Jt]]^] . Contrary to 
our dynamic treatment, their adiabatic approach models 
static disorder, which enhances the width but does not 
result in a second resonance. Hence we conclude that 
phonon-magnon interaction processes cannot explain the 
magnetic MIR and Raman anomalies. 

A severe constraint for any interpretation of the high 
frequency spectral weight is the observed increase by a 
factor of more than 2 from 4 to 300 K (triangles in inset of 
panel 2 of Fig. 2). A similar behavior can be detected in 
the temperature dependence of a(oj) of Sr 2 Cu02Cl2 @. 
We propose that these findings support the notion of a 
strong local deviation from the Neel state which is even 
more pronounced for 300 K. Whereas the broken symmc- 
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try of the antiferromagnetic state will still support long 
wavelength spin-wave excitations, the character of the 
short wavelength magnetic excitations reflects the strong 
quantum fluctuations and consequently they are insuf- 
ficiently represented by spin- waves. Both Raman and 
MIR are dominated by short wavelength magnetic ex- 
citations which makes their evaluation within spin-wave 
theory less reliable. We emphasize that this interpreta- 
tion does not contradict the good agreement of neutron 
scattering results with spin-wave theory for small mo- 
menta. Note that even in ID the inapplicability of a 
spin-wave picture to neutron data was not realized for 
many years 35 1. 

Exact diagonalization should serve to identify the 
weight of magnetic excitations. However, only clusters 
of up to V20 x V20 were investigated and produced a 
minor contribution to MIR absorption at high frequen- 
cies 36 1 . Since this cluster size is still comparable to the 
size of the considered high frequency excitations, a fi- 
nite size scaling analysis would be a serious test whether 
the Heisenberg model or a 4-spin extension || of it can 
generate the measured high frequency weight. 

In the absence of doping, the only alternative to a mag- 
netic origin of the anomaly are excitons. Perkins et al. 
jl2| suggested a d-d exciton with sidebands. However, 
this exciton should not be IR active, and there is the- 
oretical |37],|3^] and experimental |3!| evidence that its 
actual energy is a factor of two to three higher. Wang 
et al. p[ predict a CT exciton at 0.8 eV (6500 cm" 1 ) 
from fits to EELS data between 2.5 and 4 eV. We con- 
sider the Coulomb attraction necessary to pull this ex- 
citon down to 3800 cm -1 - far below the CT gap - as 
unrealistically large. Furthermore, a CT exciton should 
follow the strong redshift of the onset of CT absorption 
with increasing temperature (10500 to 9000 cm -1 , see 
Fig. 1), which is not observed. Both exciton models fail 
to describe the increase of spectral weight with temper- 
ature. Hence the exciton interpretation is an unlikely 
scenario for the considered resonance. Direct excitation 
of 2 magnons in bilayers via spin-orbit coupling was sug- 
gested [|flj to reproduce the MIR lineshape in YBC06 up 
to 4000 cm -1 . However, a spin exchange of J12 = 0.5 J 
had to be assumed, which is not anymore consistent with 
neutron scattering p5| , p6| . 

We conclude that present day understanding of mag- 
netic excitations in undoped cuprates is not sufficient to 
explain both MIR and Raman data. Only a more appro- 
priate treatment of the short wavelength excitations will 
explain the observed anomalies. 
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